Upper Bound on the region of Separable States near the Maximally Mixed State 
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A lower bound on the amount of noise that must be added to a GHZ-like entangled state to 
make it separable (also called the random robustness) is found using the transposition condition. 
The bound is applicable to arbitrary numbers of subsystems, and dimensions of Hilbert space, and is 
shown to be exact for qubits. The new bound is compared to previous such bounds on this quantity, 
and found to be stronger in all cases. It implies that increasing the number of subsystems, rather 
than increasing their Hilbert space dimension is a more effective way of increasing entanglement. 
An explicit decomposition into an ensemble of separable states, when the state is not entangled,is 
given for the case of qubits. 



I. INTRODUCTION 

A key distinguishing feature of quantum physics from 
classical physics is the prediction of a new kind of corre- 
lation between physical quantities, called entanglement. 
Quantum entanglement has been often been referred 
to as the inseparability of composite quantum systems. 
Such an entangled composite system is said to be insepa- 
rable because it cannot be prepared by manipulating each 
subsystem separately, using only measurements and op- 
erations local to one subsystem at a time. If a composite 
quantum mechanical state is specified by some density 
matrix, how can we tell if the system is entangled? 

Much work has been done studying the particular case 
of two subsystems, with each in a two dimensional Hilbert 
space (qubits). There is a good understanding of the 
entanglement for such systems and in fact a criterion, 
the partial transposition condition of Peres , indicates 
whether the subsystems are entangled. This, however, 
is a necessary and sufficient condition only when there 
are two subsystems, one with Hilbert space dimension 
2, and the other of dimension 2 or 3, as was shown by 
the Horodeckis For more complex systems, it only 
determines whether the state contains distillable entan- 
glement, however there are also some states with bound 
entanglement, which are lumped together with the sepa- 
rable states by this criterion. 

Lewenstein et. al. ^ have used the Peres condition to 
consider two subsystems, but with each subsystem now 
in a A'^ > 2 dimensional Hilbert space. Zyczkowski et. al. 

have, among other results, shown that all the mixed 
states in a sufficiently small neighbourhood of the maxi- 
mally mixed state are separable. They also gave a bound 
on the size of this neighbourhood for small composite sys- 
tems. Vidal and Tarrach ||] gave a lower bound on the 
size of this neighbourhood for arbitrary composite states, 
of any number of subsystems. Schack and Caves Q gave 
bounds for composite systems composed of many qubits 
(^ = 2). 

It has been pointed out by Braunstein et. al. [^j that 



Maximally-entangled states of the GHZ type, with noise 
added, are connected to recent proposals for NMR quan- 
tum computing. They are also relevant to fundamental 
tests of quantum mechanics using Bell inequalities. The 
separability of such noisy generalised singlet states has 
been considered by various authors recently, mostly for 
the case of qubits (subsystems with Hilbert space dimen- 
sion two). Schack and Caves using the approach of 
Braunstein et. al. have obtained an exact boundary con- 
dition for the Werner states |^ (two qubits), while Caves 
and Milburn |^ extended the approach to two q-trits 
(Hilbert space dimension three). The Horodeckis 0] 
gave exact bounds for the case of two subsystems. 

In this paper we extend the approach of Peres to 
consider such noisy generalised singlet states in the gen- 
eral case of D subsystems with each subsystem in a 
N dimensional Hilbert space. A parameter e specifies 
the amount of the pure GHZ-type maximally-entangled 
states present compared with the maximally-mixed noise 
state. We then ask and answer the following question: 
what is the maximum value of e for this state to be entan- 
gled? Another fundamental question is also considered. 
To create as entangled a state as possible, is it simply 
better to increase the dimension of the subsystems or is 
it better to increase the number of subsystems? 



II. GENERALISED WERNER STATES 

The states considered here, consist of a mixture of 
the maximally-mixed noise states, and the GHZ-type en- 
tangled states, where the number D of subsystems over 
which entanglement occurs, and the Hilbert space dimen- 
sion N of each subsystem (equal for all D subsystems) 
can take on any values iV > 1, _D > 1. The relative 
proportion of GHZ-type states is controlled by the pa- 
rameter £, where e — 1 corresponds to a pure GHZ-type 
state, while e = corresponds to the maximally mixed 
state. Explicitly, the state is given by the density opera- 
tor 
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Here \ik) represents one of N complete orthogonal ba- 
sis states for the fc"^ subsystem, while /5„ is the density 
operator for the maximally mixed state, and pe for the 
GHZ-hke state. 

The Werner state |^ , consisting of a singlet state and 
some noise, is the simplest D = N = 2 case, hence we 
call these generalised Werner states. The state p can 
be viewed as a generalised singlet state pe after it has 
emerged from a depolarising channel. 



III. SEPARABILITY 

The separability of particular cases of the state (|l|), and 
of more general states, has been considered previo usly by 
a number of authors . Firstly, it was shown [ pl]Jl2| 

that for the two-qubit {N = 2, D — 2) case, (The Werner 
state) p is separable for e < 1/3, and entangled other- 
wise, whereas for three qubits (TV = 2,D = 3), Schack 
and Caves Q found the state to be separable for e < 1/5. 
Soon after, Horodecki and Horodecki |l^ found the exact 
result for arbitrary numbers of qubits: 
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Schack and Caves also found bounds on the size of the 
separable neighbourhood around the maximally mixed 
state for totally general states of many qubits {N = 2). 
For Z? = 4 and D = 5 these are e < 1/33 and e < 1/243 
respectively, and for higher D are given by 
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What about the more general case when N and D are 
arbitrary? For what values of e is the state given by 
(|l|) separable? Vidal and Tarrach [|| gave a maximum 
bound for the random robustness R of arbitrary multi- 
component states. For the states considered here, the 
critical value of e at which the states change from being 
entangled to separable is Sc = 1/(1 + R{pe\\pn)) (in the 
notation of ref. ||]). So according to that bound. 
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We will prove in the next section (section IV) that the 
states given by (□), are always entangled if 
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That this bound is strong for qubits {N = 2) is shown in 
appendix (|^, and an explicit decomposition into prod- 
uct states is given for the case of separable p. 



IV. OUTLINE OF THE PROOF OF (6) 

Peres |^ has shown that a necessary condition for a 
state consisting of two subsystems to be separable is that 
the partial transpose of the density matrix over one of 
the subsystems, and the partial transpose over the other 
subsystem, have positive eigenvalues. However, this is a 
necessary and sufficient condition only when one of the 
subsystems has Hilbert space dimension 2 or 3, and the 
other dimension 2, as was shown by the Horodeckis 
That paper went on to give a necessary and sufficient con- 
dition for a state to be separable. Nevertheless the Peres 
condition is just what is needed for an upper bound on e 
for separable states, i.e. any states which break the con- 
dition are entangled, although some which satisfy it may 
also be, but do not have to be, entangled. It is worth 
noting that any states which satisfy the Peres condition 
but are entangled are said to contain only "bound" en- 
tanglement, as it cannot be used for teleportation, nor 
distilled by the process of entanglement distillation. 

Firstly, note that the Peres condition is easily extended 
to more than two entangled subsystems. If there are D 
subsystems, one simply chooses some group of M < D 
original subsystems to be called half-system number 1, 
and the remaining subsystems to be called half-system 
number 2. If for any such group of subsystems, an eigen- 
value of the partial transpose of p over half-system num- 
ber 1 (say) is negative, then p is entangled. Thus to use 
the Peres condition to full advantage, one must consider 
all such groups of subsystems. 

The state (|) is convenient in this respect, because it 
is unchanged under relabeling of the subsystems (evident 
by inspection). Thus the eigenvalues of pj, the partial 
transpose of p over the set a of subsystems, need only 
be looked at for D/2 (rounded down) sets of subsystems 
to extract the maximum benefit from the Peres condi- 
tion. In particular, a choice of sets of subsystems can be 
UM = {1,2,...,M} where M = 1, 2, . . . , 15/2, and an 
contains the labels of the subsystems to be considered as 
members of the Mth half-system. 

Firstly let's consider M = 1, i.e. the first subsystem's 
entanglement with the remaining D — 1 of them. The 
partially transposed density matrix is 
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Since all of the elements of p, (hence ) are finite, 
the eigenvalues of pj must also be finite. Now, exploit- 
ing the general continuity property of eigenvalues of pf , 
we conclude that at the value of e = above which the 
Peres condition indicates the state is entangled, one or 
more eigenvalues of must be zero, since they are all 
positive for e < Eq, and at least one is negative for e > So- 
i.e. for some nonzero eigenvector 
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we must have p\{eo) — 0- Expanded, this gives 
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Equation (^) can explicitly be written out as A^'^ equa- 
tions 

(1 £o)V'iii2---»D 

-I- eoA^'°~"^^i2i3^i2U • • ■ ^424DV'i2n...il = Ol (10) 

where 5ab = 1 if a = fe, otherwise. Now if one or more 
of the la '■ a — 3, . . . , D does not equal ^2, then that 
equation is satisfied only if e = 1 or ipij^i2...io = 0. The 
first case is not of interest here, as e = 1 corresponds to 
our maximally entangled GHZ like states, so we choose 

'0ii,i2---»D ~ 0. 

The rest of the equations where J2 = «3 = • • • = in, 
separate into iV(iV — l)/2 coupled sets of two equations 
of the identical form 

(1 - £o)^afc...b + eoiV""" V&a...a = 0, (11) 
(1 - eo)^Pba...a + £oN''-^^Jab...b = 0. (12) 

These have solutions if ipab...b — ^ba...a = 0, but this 
would imply \(j}) =0, which was specifically excluded in 
(^. Otherwise, these coupled two equations are only sat- 
isfied if 
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So for 1 ^ 0, at least one such coupled set of two equa- 
tions leads to the expression (p^). This is the only can- 
didate for the point where the Peres condition becomes 
satisfied. 

Now it can be easily seen that in the total-noise case 
e — 0, pj = p and all the eigenvalues of pi are 1/{N^). 



In the no- noise case {e = 1), proceeding in similar fashion 
to before, the eigenvalues A of must satisfy 
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This gives ^Hii2...iD = if for some a = 3,4, 

^2 7^ ia^ and leads to sets of two coupled equations of 

the form 

1pab...b = >^N%j}ba...a, (15) 
i'ba...a = >^Nlpab...b- (16) 

These have the solutions 'ipab...b — '4'ba...a — or 

(17) 

As before, for nonzero eigenvectors, the first cannot be 
true. So, finally, for the no-noise state, at least one eigen- 
value of p\ must be negative, and equal to 
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Thus, finally, since Eq is the only value of e where an 
eigenvalue of p\ is zero, all eigenvalues are positive for 
e = 0, and an eigenvalue is negative for £ = 1, there must 
be at least one negative eigenvalue for e > Eq given by 

m- 

When one proceeds in the same fashion when M ~ 
2, 3, . . . , Z5/2, one always gets coupled sets of two equa- 
tions identical in form to (ITTI), as can be seen by in- 
spection, so nothing new is found. Thus the result of 
section (O) is indicated. 



V. COMPARISON TO KNOWN BOUNDS 

For the qubit case (A^ = 2), as seen in figure (|l|), equa- 
tion (H) gives an exact bound on the values of e that di- 
vide separable from entangled states of the form (^. One 
sees that the upper bound (|^) derived from the work of 
Vidal and Tarrach comes very close to the exact value 
for the qubit case. This exact bound is also greater than 
those lower bounds previously found by Schack and Caves 
Q, as expected. 

For the two-subsystem [D = 2) case, the bound (H) 
agrees with the exact one found by the Horodeckis |10( |. 

For other values of D and iV, the results on random 
robustness by Vidal and Tarrach, lead to an upper bound 
which is considerably weaker than the upper bound (|^) 
given by the partial transposition condition. The upper 
bound found here, actually gives a stronger bound on 
the random robustness of entanglement of states given 

by (0). 
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R{Pe\\Pn) < N''-^ (19) 

It is interesting to note that in the border cases when 
TV = 2 or D = 2, the bound (|) is in fact an exact bound. 
This is despite the fact that the Peres condition does not 
necessarily give a strong bound for such states. This leads 
one to the tentative conjecture that for noisy GHZ-type 
states of the form (|l|) , the Peres condition may in general 
give a strong upper bound. 

Looking at figure (||) , one sees that as the Hilbert space 
dimension of the subsystems increases, the upper bound 
on £ rapidly decreases, indicating that the entanglement 
becomes stronger. 

As the bound (||) is completely general in N and D, 
it does provide some answers to the question of what 
raises entanglement more: creating more entangled sub- 
systems, or increasing their dimension? Since the bound 
on e decreases exponentially with D, but only polynomi- 
ally with TV, one concludes that increasing the number 
of subsystems is a more effective way of increasing the 
entanglement. 



APPENDIX A: PROOF THAT (6) IS EXACT FOR 
QUBITS 

We wish to show that the bound (||) is exact for qubits 
{N = 2), and to find the product states which combine 
to give the state when it is separable. 

Starting off similarly to Schack and Caves , the D- 
subsystem generalisation of the Werner state ((^ with 
N = 2) can be written in terms of Pauli matrices: 

Ks)^^{{l~e)I^^+'-E}, (Ala) 

where 

E= {I + as)'^^ + {I - <J3)^D 

+ {<Jl + iT2)®^ + (fTi - 1^2^^, (Alb) 

and ai are the Pauli matrices, / is the two-dimensional 
identity matrix, and for conciseness, the following nota- 
tion is used: 



VI. CONCLUSION 

The results presented here based on the Peres condi- 
tion provide a lower bound on the parameters e in (0) 
above which the generalised Werner states are always 
entangled. Furthermore, it gives an exact bound on this 
parameter for the case of many qubits. An explicit sim- 
ple expression is derived that depends on D, the number 
of subsystems over which entanglement occurs, and N 
the Hilbert space dimension of each subsystem. Apart 
from the few cases {D = 2;N — 2 and D = 3) where this 
bound was known exactly previously, the new bounds are 
stronger than previously known ones. 

This work also sheds light on the question of whether 
to increase quantum entanglement in a system, is it bet- 
ter to create more entangled subsystems, or to increase 
the dimension of the existing subsystem. As the bound 
on £ decreases exponentially with TV, but only polyno- 
mially with D, increasing the number of subsystems is a 
much more effective way of increasing entanglement. 

To conclude, the Peres partial transposition condition 
has provided a good upper bound for determining the 
separability of a generalised N, D Werner state. How- 
ever for other systems it is known that this transposition 
condition fails to give strong results, thus when this con- 
dition is useful, and when not, remains an interesting 
question. 
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(A)^^ = (A) • • • ® (A) D times. (A2) 

To show that (^) is actually a strong bound, is suffices 
to find an expansion of p in terms of a positive sum of 
direct tensor products of density matrices for this value 
of 



By analogy with the results of Schack and Caves [D, 
it has been guessed that an expansion of p(£c) in terms 
of the density matrices 

P±, = i(/±a,) (1 = 1,2,3), (A4) 

is given by the following: 

(A5) 

Here the sum is over all permutations of D indices 
xi, X2, ■ ■ ■ , xd satisfying the following conditions: 

1. e {1,-1,2,-2}. 

2. The number of Xi G {2, —2} is even (or zero). 

3. If the number of Xi G {2, —2} is a multiple of four 
(or is zero), then the munber of .t^ G { — 1, —2} is 
even (or zero). 

4. If the number of Xi G {2,-2} is not a multiple of 
four, then the number of Xi G { — 1, —2} is odd. 
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The proof of this proceeds by starting with the above 
guess, and showing that it is a correct one. As it turns 
out, the hard work is in actually writing down the guess 
mathematically, so let us begin with this. 

To begin with, note that 

ro = (Pi+P2+P-l+P-2)®^ 



P:., ® (A6) 



where the sum is over all permutations of D indices 
xi,X2, ■ ■ ■ ,xn & {1, 2, —1, —2}. Also note that 



So = (Pi + P2 - P-1 - P-2 



(A7) 



will give a sum over all these permutations, except 
that all terms which have an odd number of indices in 
{—1,-2} will be subtracted rather than added like in Tq. 
Thus one can see that 



(A8a) 



will give a sum over permutations of D indices, like Tq, 
except that only terms where the number of indices in 
{ — 1, —2} is even will be included. Similarly 



^0 — 2^'^ ^ '^"1 



(A8b) 



will give a sum over only those terms in which the number 
of indices in { — 1,-2} is odd. 

Now consider some more similar expressions. 



'ri = (Pi-P2 + P-i-P-2)^''. 



(A9) 



71 will give a sum over all index permutations, except 
that all terms which have an odd number of indices in 
{2, —2} will be subtracted rather than added like for Tq. 



72 = {Pi + 1P2 + P-i + iP-'. 



(AlO) 



T2 will give a similar sum over all permutations, but 
terms in which the number of indices in {2, —2} is a mul- 
tiple of four (or is zero) will be added, terms in which 
this number is even, but not a multiple of four, will be 
subtracted, terms in which this number is one more than 
a multiple of four will be added and multiplied by i, and 
terms for which this number is one less than a multiple 
of four will be subtracted and multiplied by i. 



% = (Pi -1P2 + P-1- iP-: 



(All) 



73 is the complex conjugate of 72. It can be seen (after 
a little thought) that -j [7o + 7i + 72 + 73] will give a sum 
over only those terms in which the number of indices 
in {2, —2} is a multiple of four (or is zero). Similarly, 
I; [7o + Ti — 72 — 73] will give a sum over only those terms 
in which the number of indices in {2, —2} is even, but 
not a multiple of four. 



Analogously to equation (A7) define 



Si = (Pi - P2 - P-i + P-2)®^, (A12a) 

52 - (Pi + iP2 - P-i - *P-2)®^, (A12b) 

53 = (Pi - *P2 - P-l + lP^2f, (Al2c) 

and one can define expressions fo r TZ f and TZ° for i = 
0, 1, 2, 3 analogously to equations (A8). So following the 
same reasoning as previously, the sum of all terms in 
which the number of indices in {2, —2} is a multiple of 
four, and the number of indices in { — 1,-2} is even is 
given by 



(A13) 



And thus finally, the sum in the second term of the 
guess pg ( equation ( [A^ ) ) can be written 



1, 



(A14) 



[To + Ti + 52 



So the guess that has been made (i.e. equation (|A5|) ) can 
be rewritten 



P 



2N+ 



J [To +T1+S2+ S3] 



(A15) 



which is explicitly (via the expressions (A6), ([A7|), (|AS|), 
( A12| )) a positive sum of direct tensor products of density 
matrices, and thus is separable. The only question that 
remains is whether pg = p(£c) ? 

This is the easy part. It is seen using the expression 



(A4) that 



To 
Ti 
S2 
S3 

^3 



0, 

(ai+*'T2)«^, 
(ai-*'T2)«^, 

2-^(7 + fT3)^^, 



p^iv ^ 2-'' {I - a3f- 



so 



2N+I 



(Al6a) 
(A16b) 
(Al6c) 
(A16d) 
(A16e) 
(A16f) 



(A17) 



which only seemingly differs from equation (Al) by the 
first term, but using the expression for Ec (A3), one finds 
that these first terms are equal also. 



1 



2 ' 



(A18) 



so Pg = p, the guess was correct, and thus the bound Sc 
is strong. 
Q£V 
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FIG. 1. Bounds on the value of e for which the states (|l|) 
become separable or entangled, for the qubit (A*' — 2) and 
q-trit (A'' = 3) cases, shown on the same plot. When e is 
above the upper bounds or the exact bound, p is entangled, 
and when e is below the lower bound or the exact bound, p 
is always separable. When e is below an upper bound, p may 
be separable or bound entangled. 




FIG. 2. Upper bounds on e for separable states of the 
form (jl|) on a logarithmic scale. Variation with subsystem 
Hilbert-space dimension A*' is shown. 
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